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Enumerating nonisomorphic orthogonal arrays is an important, 
yet very difficult, problem. Although orthogonal arrays with a spec- 
ified set of parameters have been enumerated in a number of cases, 
general results are extremely rare. In this paper, we provide a com- 
plete solution to enumerating nonisomorphic two-level orthogonal ar- 
rays of strength d with d + 2 constraints for any d and any run size 
n = A2 d . Our results not only give the number of nonisomorphic or- 
thogonal arrays for given d and n, but also provide a systematic way 
of explicitly constructing these arrays. Our approach to the problem 
is to make use of the recently developed theory of J-characteristics 
for fractional factorial designs. Besides the general theoretical results, 
the paper presents some results from applications of the theory to or- 
thogonal arrays of strength two, three and four. 

1. Introduction. An orthogonal array of size n with m constraints, s 
levels and strength d > 2 is an n x m matrix with entries from a set of s 
levels, usually taken as 0, 1, . . . , s — 1, such that for every n x d submatrix, 
each of the s d level combinations occurs the same number A of times. Such an 
array is denoted by OA(n, m, s, d). The five parameters n,m, s,d, A cannot 
vary independently, as the relation n = Xs d must hold. This is why A, called 
the index of the array, is suppressed in the notation OA(n, m, s, d). 

The importance of orthogonal arrays cannot be overstated. They are in- 
timately related to many other combinatorial objects such as Hadamard 
matrices, orthogonal Latin squares and error-correcting codes. They are di- 
rectly useful as fractional factorial designs in factorial experiments arising 
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from scientific and technological investigations. The most fundamental ques- 
tions concerning orthogonal arrays are when they exist and how to construct 
them if they do exist. Indeed, these are the two main themes of [13]. Con- 
struction and optimality of orthogonal arrays as fractional factorial designs 
are discussed in [9] . A comprehensive discussion on the use of orthogonal ar- 
rays as fractional factorial designs is given by Wu and Hamada [23]. When 
more than one orthogonal array is available for a given factorial experiment, 
the experimenter can use a suitable criterion to choose among the candi- 
date arrays. The most commonly used criteria are those of resolution and 
aberration, along with their generalizations [4, 10, 16, 20, 26] and [24]. 

This paper considers the enumeration of orthogonal arrays, an important, 
yet very difficult, problem. For given parameters, finding the complete set 
of orthogonal arrays is of theoretical importance in its own right. Knowing 
the complete set of orthogonal arrays allows us to randomly select one to 
use from all arrays, an important consideration in designing an experiment. 
In addition, such knowledge of all orthogonal arrays often provides useful 
experimental data for theorists. Furthermore, having the complete set of 
arrays available may lead to an unexpected discovery of orthogonal arrays 
with other useful statistical properties besides those directly offered by the 
strength of these arrays. Enumerating orthogonal arrays can be greatly sim- 
plified by considering only nonisomorphic orthogonal arrays, as all arrays 
can easily be generated once all nonisomorphic arrays are available. Two 
orthogonal arrays are said to be isomorphic if one can be obtained from the 
other by a sequence of operations involving permuting rows, columns and 
levels, and are said to be nonisomorphic otherwise. 

Despite the importance of the problem, literature on the subject is con- 
siderably scarce as compared to that on the existence and construction 
problems, due to the obvious reason that the enumeration problem is ex- 
tremely difficult. Here, we survey some main findings. Lam and Tonchev 
[14] enumerated all 04(27, 13, 3, 2) 's. Hedayat, Seiden and Stufken [12] com- 
pletely classified all 04(54, 5,3, 3) 's. Saturated two-level orthogonal arrays 
04(4A,4A - 1,2,2) for A = 1,2,..., 6 are enumerated in [27, 28] and [13], 
Theorem 7.37. Sun, Li and Ye [21] provided a complete enumeration of 
04(16, m, 2, 2)'s for all 3 < m < 15 and 04(20, m, 2, 2)'s for all 3 < m < 19. 
Recently, Bulutoglu and Margot [2] completely enumerated 04(80,6,2, 4) 's, 
04(96,7, 2, 4)'s, 04(112, 6, 2, 4)'s and 04(144, 8, 2, 4)'s. Chen, Sun and Wu 
[5] and Xu [25] considered enumeration of linear orthogonal arrays (regu- 
lar fractional factorial designs) for small run size n. Yumiba, Hyodo and 
Yamamoto [30] classified all 04(24, 6, 2, 2)'s that can be obtained from sat- 
urated 04(24, 23, 2, 2)'s. The above enumeration results are all for specific 
sets of parameters. The most general result we know of is that of Seiden and 
Zemach [17], which gives a complete enumeration of two-level orthogonal 
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arrays of strength d with d+1 constraints. Two other results of some gener- 
ality were reported by Fujii, Namikawa and Yamamoto [11], who completely 
classified OA(2 d+1 ,d + 2,2,d)'s and OA(2 d+2 , d + 3, 2, cf)'s. Orthogonal ar- 
rays are special cases of balanced arrays. Srivastava [19] and Shirakura [18] 
presented the numbers of nonisomorphic balanced arrays of strength d for 
m = d+l,d + 2 and d + 3 factors. However, their definition of isomorphism 
is different, in that only row and column permutations are used. In fact, it 
would not make sense to consider symbol permutations in defining isomor- 
phic balanced arrays, as permuting the symbols of a column in a balanced 
array may result in an unbalanced array. 

In this paper, we provide a complete solution to the problem of enu- 
merating two-level orthogonal arrays of strength d with d + 2 constraints. 
The results are general in that they are for any strength d and any in- 
dex A. Our approach is to make use of the recently developed theory of 
J-characteristics for fractional factorial designs [22]. The remainder of the 
article is structured as follows. Section 2 reviews J-characteristics, discusses 
their application to orthogonal arrays and characterizes nonisomorphic ar- 
rays using J-characteristics. Section 3 presents our main results through 
four theorems. Applications to orthogonal arrays of strength two, three and 
four are also discussed here. Section 4 concludes the paper with a summary 
of our findings and a discussion of possible future work. 

2. Enumeration method. Our method of enumerating two-level orthog- 
onal arrays is based on the recently developed theory of J-characteristics 
for fractional factorial designs. Section 2.1 briefly reviews fractional facto- 
rial designs and their J-characteristics, and presents a key identity linking 
a fractional factorial with its J-characteristics. Section 2.2 specializes J- 
characteristics for orthogonal arrays of strength d with d + 2 constraints. 
Section 2.3 considers how to characterize nonisomorphic orthogonal arrays 
of strength d with d + 2 constraints using J-characteristics, paving the way 
for the main results to come in Section 3. 

2.1. Fractional factorial designs and their J-characteristics. For any s C 
Z m = {!)•••) m }> define a row vector 



where r s j = — 1 if j £ s and r s j = +1 otherwise. A complete 2 m factorial can 
then be represented by the 2 m x m matrix 



(1) 





C = (hi,...,h m ). 
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The common component of hj and r s is r s j, as in (1). By the sth row of 
C, we mean r s and, similarly, by the sth component of hj, we mean r s j. 
For any t C Z m , the Hadamard product hj of the columns hj with j € t is 
a column vector of length 2 m , the sth component of which is Ylj^t r sj- Let 
H be the Hadamard matrix of order 2 m formed by collecting in Yates order 
all of the Hadamard products of the columns of C, that is, 

(3) H = (h^, hi, h 2 , b.12, h 3 , hi 3 , h 23 , hi23, h 4 , h i4 , . . .). 

Deleting h^, from H gives a saturated factorial in m independent columns. 

A two-level fractional factorial design of n runs for m factors is annXm 
matrix D = (dij) with dij = ±1, where each row of D corresponds to a run 
and each column to a factor. A row permutation of D gives the same design. 
Let N s denote the number of times that a run r s in (1) occurs in D. This very 
N s is called the indicator function by some authors; the reader is referred 
to [29] for a more extensive discussion. Design D can then be equivalently 
described by the vector of length 2 m , 

(4) N = (AV, N U N 2 ,N 12 , N 3 ,N 13 , N 23 ,N 123 ,N 4 , N u , . . .) T , 

where, for example, N\ 2 is shorthand notation for iV{ 12 }- As a matrix, D is 
determined by N up to a row permutation and, as a design, D is completely 
determined by N. Conversely, any vector N of length 2 m with nonnegative 
integers as its components defines a design with n = J2scz m N s runs. 
For any t C Z m , let 

n 

(5) J t = Y^l[d ij = h st N s , 

i=l j£t s<ZZ m 

where H = (h s t) is the Hadamard matrix in (3) and h s t denotes the sth 
entry of h^, the tth column of H. The Jf values defined in (5) for all subsets 
t of Z m are called the J -characteristics of design D. Let 

(6) J = (J(j>, Ji,.h, J12, <h, Ji3, J23, J123, J4, Ju, ■ ■ -) T • 

The set of J-characteristics as in (6) completely determines a general facto- 
rial, just as the defining relation does for a regular factorial. This is given in 
the following lemma. 

Lemma 1. N = 2~ m HJ, that is, 

(7) N s = 2~ m Y, hstJt, 

tC-Zfn 

where H = (h s t) is as given in (3). 
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This result is given in [22] and plays a key role in this paper. As discussed 
in [22], J-characteristics directly capture the projection properties of a fac- 
torial design and, because of this, they lead to successful generalizations of 
minimum aberration from regular factorials to general factorials [7]. Further 
work on generalized minimum aberration can be found in [26] and [3]. Fi- 
nally, we note that, mathematically, J is simply the Hadamard transform of 
N. 

2.2. J-characteristics for orthogonal arrays. For two-level orthogonal ar- 
rays, it is convenient to use ±1 to denote the two levels. A general two-level 
fractional factorial design D = (dij) nxm , where dij = ±1 as introduced in 
Section 2.1, may or may not be an orthogonal array of strength d>2, but, 
if it is, then its J-characteristics have some special properties. The next three 
lemmas are devoted to these properties of J-characteristics of an orthogonal 
array. 

Lemma 2. Design D is an orthogonal array of strength d if and only if 
Jt = for all t C Z m with \t\ < d, where t ^ (ft. 

Lemma 2 can easily be verified using (5) and (7) and was given earlier in 
[7]- 

Lemma 3. Suppose that D is an OA(n, m, 2, d), where n = \2 d and m > 
d + 2 . Then: 

(i) for any t C Z m , it holds that Jt = [it2 d for some integer fitj 

(ii) if A is even, then fi t is even; 

(hi) if A is odd and d is even, then \it is odd for |t| = d + l,d + 2; 
(iv) if A is odd and d is odd, then fit is odd for \t\ = d + 1 and even for 
\t\ =d + 2. 

Proof. From (1), we see that = (1, . . . , 1), implying that h^t = 1 for 
any t. Taking s = (ft in (7) gives 

(8) iV = 2~ m £j 4 . 

For given w C Z m , consider the projection design consisting of the 
columns of D in w. Now applying (8) to this projection design D w , we 
obtain = 2 - ' 1 "' J2tcw Jt- From Lemma 2, we have Jt = for all t with 
1 < \t\ < d. Further, note that J^ = n = X2 d . Therefore, 

(9) J w = 2^% - \2 d - J t- 

tcw,d+l<\t\<\w\-l 
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Based on (9), Lemma 3 can easily be established. □ 

For the special case d = 2, results similar to those in Lemma 3 were pro- 
vided in [8]. Lemmas 2 and 3 discuss the possible values of Jt for a given t. 
The following result states that the values of Jt for two different subsets t\ 
and t2 are, in fact, related, provided that the two subsets are close to each 
other. 

Lemma 4. Let D be an 0A(n,m,2,d). We must have 

\JtA + \ Jt 2 \ < n , 

provided that 1 < |tiAt2| < d, where t\At2 = (t\ \ t^) U (£2 \ t\). 

Proof. Let dj be the jth column of D. Let dt 1 be the Hadamard 
product of the columns dj for j Et\. Define dt 2 similarly. Consider the design 
D* given by the two column vectors d tl and d t2 , that is, D* = (d tl ,d t2 ). 
Obviously, the J-characteristics for column vectors d^ and dt 2 are Jt x and 
Jt 2 , respectively, and the J-characteristic for the two combined is Jt 1 At 2 - 
Because 1 < |iiAt2| < d and D is an orthogonal array of strength d, we 
must have Jt 1 At 2 = 0- We apply Lemma 1 to design D*. Note that the 
corresponding Hadamard matrix H in (3) is now given by 

H= [(1,1,1,1) T ,(1,-1,1,-1) T ,(1, 1,-1, -1) T ,(1, -1,-1, 1) T ] T 

Noting that N s in Lemma 1 is nonnegative, we obtain 

n + J tl + Jt 2 > 0, n— J tl + Jt 2 > 0, 
n + Jtj — Jt 2 > 0, n — Jt x — Jt 2 > 0. 
This shows that \Jt x \ + | Jt 2 \ < n. □ 

The seemingly unremarkable result in Lemma 4 is actually very powerful. 
For example, a key result in [6], Lemma 2.2, is a direct consequence of 
our Lemma 4. In what follows, we document an interesting application of 
Lemma 4 to orthogonal arrays of strength 2, the most useful situation in 
the statistical design of experiments. 

Corollary 1. Let D be an OA(AX,m,2,2) with odd A. Then we must 
have \ Jt\ < 4A — 4 for all t with 1 < \t\ < m — 1. For t = Z m , the possibility 
that \Jt \ = 4A can occur only when the two conditions simultaneously hold 
that (i) m = 3 (mod 4) and (ii) D cannot be embedded into an 0^4(4A,m + 
1,2,2). 
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We omit the proof, which involves using a result in [8], Proposition 1, and 
is otherwise straightforward. If \Jt \ = 4A, then the set of columns dj with 
j € t is said to form a defining word, which results in full effect aliasing. 
Corollary 1 effectively states that when the run size n = 4A is not a multiple 
of 8, then, except in rare cases, an orthogonal array of strength 2 does not 
have a defining word. 

2.3. Characterizing nonisomorphic OA(n, d+2, 2, d)'s using J -characteris- 
tics. As discussed in Section 2.1, a fractional factorial design D is com- 
pletely defined, up to a row permutation, by the vector N in (4). Lemma 1 
asserts that N is uniquely determined by J in (6), and vice versa. Therefore, 
up to a row permutation, a factorial design D is completely characterized 
by its J-characteristics. When design D is an OA(n, d + 2, 2, d), Lemma 2 
asserts that = n, and Jt = for all nonempty t with \t \ <d. So, the vector 
J in (6) can be equivalently described by 

(10) J = (J tl ,...,J Wl ), 

where tj = Z m \ {m + 1 — j} for j = 1, . . . , m and t m +i = Z m . For convenience, 
we use the same notation J for this shortened version of the original vector 
J. No confusion will arise, as the new vector J in (10) is always used in the 
context of OA(n, d + 2, 2, d)'s. We call the vector J in (10) the J-vector of 
the OA(n, d + 2, 2, d) under consideration. Two orthogonal arrays are said 
to be statistically equivalent if one can be obtained from the other by a 
row permutation [13]. Thus, two statistically equivalent OA(n, d + 2, 2, d)'s 
have the same J-vector and two OA(n, d + 2, 2, d)'s that are not statistically 
equivalent have different J-vectors. The preceding discussion leads to the 
following conclusion. 

Lemma 5. There is a one-to-one correspondence between the classes of 
statistically equivalent OA(n, d + 2,2, cQ's and their J-vectors. 

It is tempting to use the J-vector to identify nonisomorphic arrays, but 
in its present form, the J-vector does not serve this purpose. This is be- 
cause two isomorphic arrays will have two different J-vectors if they are not 
statistically equivalent; for a given array, both permuting its columns and 
switching the signs of its columns affect its J-vector. Lemma 5 has taken care 
of isomorphism due to row-permuting. We now consider isomorphism due to 
column-permuting and sign-switching. For a given array, there are m! ways 
of permuting its columns. Let 5j be the indicator of whether the jth col- 
umn is sign-switched, with 8j = — 1 denoting that there is sign-switching and 
5j = +1 denoting that there is no sign-switching. A vector S — (^i, ■ ■ ■ ,^m) 
then gives one way of sign-switching the columns of the array. Thus, the total 
number of ways of sign-switching the columns is 2 m . Starting from a given 
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array, column-permuting and sign-switching together generate m\2 m arrays 
in all, some of which may be statistically equivalent. From the J-vectors of 
these m!2 m arrays, we will choose a unique J-vector, and therefore a unique 
array, to represent the whole class of isomorphic arrays. 

We first present a lemma which will be needed in what follows and also in 
Section 3. Consider the fractional factorial design B obtained by collecting 
columns bi, . . . , b m+ i from the Hadamard matrix H in (3), that is, 

(11) B = (bi,...,b m+ i), 

where bj = h t . with tj = Z m \ {m + 1 — j} for j = 1, . . . ,m and i m +i = Z m . 
Recall that a defining word in a fractional factorial design is a set of columns 
whose Hadamard product is I, a column of all +l's. 

Lemma 6. Design B has exactly one defining word, given by b± ■ ■ ■ b m+ i = 
I when m is even and h\ ■ ■ ■ b m = I when m is odd. Therefore, for even m, 
deleting any column from B gives a full factorial design and for odd m, 
deleting any column except for h m+ i also gives a full factorial. 

The proof is simple and is thus omitted. 

We now return to orthogonal arrays. Let D = (di, . . . , d m ) be an OA(n, m, 2, 
where m = d + 2 and dj denotes the jth column of D. As discussed ear- 
lier, a vector 6 = (S±, . . . ,S m ) with Sj = ±1 gives one way of switching the 
signs of the columns. Sign-switching D by S = (<5i, . . . , S m ) results in a new 
array D s = (^ldi, ■ • ■ ,S m d m ). Obviously, we have J tj (Ds) = ^^( D ) for 
j = 1, . . . ,m + 1, where 5t j = Ilietj Since the 2 m vectors 5 = (Si, . . . , 6 m ) 
collectively form a full factorial, the resulting 2 m vectors (Jtj , . . . , 6t m+1 ) are 
exactly the 2 m runs of design B given in (11). 

Now, suppose m is even. The situation involving odd m will be dealt with 
later. 

Case (i): If |Jt m+1 | < mini<j< m | J tj |, select a S = (5i, . . . ,5 m ) such that 
Stj = —1 if Jtj > and 5^ = +1 if Ju < 0, for j = 1, . . . , m. This is possible 
because for even m, the 2 m vectors (5^ , . . . , 5t m ) form a full factorial ac- 
cording to Lemma 6. Such a choice of S = (6±, . . . , S m ) leads to an array, the 
J-vector of which satisfies 

(12) Jt,<0 for j = l,...,m. 

Clearly, permuting the columns of an array results in a permutation of 
Jt x , • ■ • , Jt m ■ It is a l so eas Y to see that any permutation of J% x , . . . , Jt m cor- 
responds to a permutation of the columns. For the array satisfying (12), we 
can permute its columns to obtain an array whose J-vector satisfies 



(13) 



J tl < ■ ■ ■ < Jt m < o. 
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The J-vector resulting from the above procedure is unique, as is its corre- 
sponding array. We note that although there may be more than one 5-vector 
or more than one permutation of columns that can lead to a J-vector sat- 
isfying (13), this J-vector is always unique. For example, if = 0, we can 
also use a 5-vector having S^ = —1 instead of St 1 = +1, which will not have 
any effect on the resulting J-vector since both Jt x and Jt m+1 are still zero. 

Case (ii): If |J tm+1 1 > min^x™ \J t . \ = \ J t .,\, select a S = (Si, ...,6 m ) such 
that S t . = -1 if J t . > and S t . = +1 if J tj < for j = 1, . . . , f - l,f + 
1, . . . , m + 1. This is possible because for even m, the 2 m vectors (^ 1 , . . . , St ., 
St., j , • • • , St m+ i ) also form a full factorial according to Lemma 6. This choice 
o/ 6 gives an array whose J-vector satisfies 

J t m+ i < and J t . < for j = 1, l,f + 1, ... , m. 

Now, permute J tl , . . . , Jt m by first moving J t ., to the end of the sequence 
and then arranging the others in ascending order. We now obtain an array 
whose J-vector satisfies 

(14) J il <---<Jt m _ 1 <0 and J Wl <0. 

Once again, the resulting J-vector is unique, as is the corresponding array. 
It is possible that more than one j' satisfies mini<j< m \Jt j \ = \Jt -,\- We can 
easily check that it has no effect on the final outcome. 

The above results can be summarized in the following proposition. 

Proposition 1. When m is even, every class of isomorphic orthogonal 
arrays of strength d with m = d+2 constraints contains a unique array whose 
J-vector satisfies either the condition in (15) or that in (16) as given by 

(15) Jt 1 <---<Jt m <-\Jt m+1 \, 

(16) Jt 1 <-<J tm - 1 <-\Jt m \, Jt m+1 < -\JtJ- 

Proposition 1 is immediate as (15) is (13) under case (i) and (16) is (14) 
under case (ii). 

The situation for odd m is considerably simpler. Let j' be such that 
\Jt.,\ = mini<j< m , | J t \. We select a S= (Si, . . . ,S m ) such that 5 tj . = — 1 if 
J t -> and S h = +1 if J tj < for j = 1, . . . , j' - l,f + 1, . . . ,m + 1. This is 
possible because for odd m, the 2 m vectors (S^ , . . . , St ., _ 1 , St ., , ■ ■ ■ , St m+1 ) 
form a full factorial according to Lemma 6. This choice of d gives an array 
whose J-vector satisfies 

Jt m+1 < and J tj < for j = 1, . . . , j' - l,j' + 1, . . . , m. 

Now, permute J tl , . . . , Jt m by first moving J t ., to the end of the sequence 
and then arranging the others in ascending order. We now obtain an array, 
the J-vector of which satisfies 

J tl <---< Jt m -!<0 and J Wl <0. 
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Using the same arguments as before, we conclude that the resulting J-vector 
is unique. 

Proposition 2. When m is odd, every class of isomorphic orthogonal 
arrays of strength d and with m = d + 2 constraints contains a unique array 
whose J -vector satisfies 

(17) J h < ■ ■ ■ < Jt^ < -\JtJ, Jt m+1 <0. 

To enumerate nonisomorphic OA(n, d + 2, 2, d)'s, Propositions 1 and 2 
conclude that we only need to consider OA(n, d + 2, 2, d)'s whose J-vectors 
satisfy (15) or (16) for even m and (17) for odd m. For convenience, a 
J-vector is called a J*-vector if it satisfies (15) or (16) for even m and 

(17) for odd m. Since a J*-vector uniquely determines a class of isomor- 
phic OA{n, d + 2, 2, d)'s, two (A4(n, d + 2,2, d)'s are nonisomorphic if their 
J*-vectors are different. We therefore establish a one-to-one correspondence 
between nonisomorphic OA(n, d + 2, 2, <i)'s and their J*-vectors. Enumerat- 
ing orthogonal arrays of strength d with m = <i + 2 constraints is equivalent 
to enumerating their J*-vectors. 

3. Enumeration results. Our results on enumerating OA(X2 d , d+2, 2, d)'s 
are presented in Sections 3.1, 3.2, 3.3 and 3.4, respectively, dealing with the 
four cases (i) even d and odd A, (ii) even d and even A, (hi) odd d and odd 
A and (iv) odd d and even A. Since these four sections are structurally sim- 
ilar, we shall only preview Section 3.1 here. The main result is contained in 
Theorem 1, which establishes a one-to-one correspondence between the set 
of nonisomorphic OA(X2 d , d + 2, 2, d)'s with even d and odd A and the set of 
solutions to equation (18) with its unknowns satisfying certain conditions. 
Example 1 and Corollary 2 give a first taste of the power of the theorem. 
Through Lemmas 7 and 8, we then provide a systematic method for explic- 
itly obtaining all solutions to equation (18). Finally, some results from the 
application of Theorem 1 and Lemmas 7 and 8 are presented in Tables 1 
and 2. 

3.1. Enumerating OA(X2 d , d + 2, 2, d)" "s for even d and odd A. Recall that 
m = d + 2. Consider the equation 

(18) A + uiH \-u m+ i=4k, 

where k and the Uj's are the unknowns, k > is an integer and the Uj's are 
odd integers with \uj | < A — 2 for all j = 1 , . . . , m + 1 . Further, the Uj 's must 
satisfy one of the conditions 

(19) ui<---<u m <-\u m+1 \, 

(20) ui < ■ ■ ■ <u m -i < -\u m \, u m+ i < —\Ufn\ — 2. 



ENUMERATING ORTHOGONAL ARRAYS 



11 



It is clear that (19) and (20) are mutually exclusive. Note that both the 
solutions to (18) under (19) and those to (18) under (20) are considered in 
the following theorem. 

Theorem 1. Suppose that d is even and X is odd. We then have that: 

(i) each solution (u\ , . . . , u m +i , k) to equation (18) under either (19) or 

(20) determines an OA(X2 d , d + 2, 2, d) with J* -vector given by Jt j =2 d Uj 
for j = l,...,m + l; 

(ii) the complete set of nonisomorphic OA(X2 d ,d + 2,2,d) , s is given by 
collecting the arrays obtained in (\) over all the solutions to equation (18). 

Proof, (i) Let (u\, . . . ,u m +x, k) be a solution to (18) with all of the 
accompanying conditions satisfied. Consider 

(21) J* = (J tl ,...,J tm+1 ), where J t . =2 d Uj for j = l,...,m + l. 

We will show that the vector J* = ( J tl , Jt m+1 ) in (21) is indeed the J*- 
vector of an OA(X2 d ,d+2,2,d). Obviously, (19) implies (15) and (20) implies 
(16). Prom Lemma 1, it is clear that our goal is achieved if we can show that 

/ m+l \ 

N s = 2- m [n+J2h stj J tj j 

is a nonnegative integer for all s C Z m . Noting that m = d + 2 and n = X2 d , 
we have 

(m+l \ / m+l \ 

A + ]T h stjUj /4=(A+£u j /4+£ [h st . - l]u 3 /A. 
j=i J \ j=i J j=i 

Since A + Y^j^i u j = ^ an< ^ h-stj = il 5 we obtain 

N s = k- uj/2 = k- A/2, 

j:h st] =-l 

where A = J2j-h st =-i u j- Since k > is an integer, if we can show that 
A is a nonpositive even integer, then N s must be a nonnegative integer. 
From Lemma 6, we see that the number of j's with h s t i = —1 has to be even 
because m is even and (h s t 1 , . . . , h s t m+1 ) is a run of design B in (11). We now 
obtain that A is even, since all Uj's are odd and so A = J2j ■. h st .=-1 u j * s a 
sum of an even number of odd numbers. For s = (ft, we have N s = k > 0. For 
a nonempty s, the number of j's with h s t j = —1 is a positive even integer, 
implying that J2j : h st =-i u j ^ s a sum °^ a ^ l eas t two Uj's. From (19) or 
(20), we see that the only Uj that can possibly be positive has the smallest 
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absolute value among all Uj's. [This possibly positive Uj is u m+ i if (19) holds 
and is u m if (20) holds.] We can now conclude that A < 0. 

(ii) Part (i) of Theorem 1 states that a solution (ui, . . . ,u m+ i, k) to (18) 
determines an OA(X2 d , d+2, 2, d) with J*-vector given by Jtj = 2 d Uj. Clearly, 
two different solutions to (18) give two different J*-vectors, representing two 
nonisomorphic OA(X2 d , d + 2, 2, <i)'s from Proposition 1. Part (ii) of Theo- 
rem 1 is established if we can show that every OA(X2 d , d + 2, 2, d) can be 
obtained from a solution to (18). Now, suppose an OA(X2 d , d + 2, 2, d) has a 
J*-vector given by 3* = (J tl ,..., Jt m+1 ). Let uj = Jtj1 d for j = 1, . . . , m + 1. 
From Lemma 3(iii), we know that the u^s are odd integers. Lemma 4 fur- 
ther implies that \uj \ + \v,j>\ < X for j ^ j' . Since Uj, uy and A are all odd, 
we must have \uj \ < A — 2 for all j = 1, . . . , m + 1. Now, consider 

/ m+1 \ 

^ = 2" m + J t ,J = (A + Ul + ■ ■ ■ + u m+1 )/4. 

We then see that (u\, . . . ,u m+ i, k), where k = N^, is a solution to equa- 
tion (18). The condition in (19) or (20) is met due to Proposition 1. □ 

For even d and odd A, Theorem 1 establishes a one-to-one correspondence 
between the set of nonisomorphic OA(X2 d , d + 2, 2, d)'s and the set of solu- 
tions to equation (18) subject to the given conditions. From the above proof, 
we see that N s > Nj, = k > 0, so the k in a solution (ui, . . . , u m +i, k) to (18) 
has a simple interpretation: the corresponding array contains k copies of a 
full factorial, with k being the largest integer for which such a statement 
can be made. 

Example 1. First, consider the case A = 3, d = 2, n = 12 and m = 4. 
Equation (18) has a unique solution under (19), given by (ui, U2, «3, «4, «5, k) = 
(— 1, -1, -1, -1, 1,0), and has no solution under (20), so the 0^(12,4,2,2) is 
unique up to isomorphism. For the case A = 5,d = 2,n = 20 and m = 4, equa- 
tion (18) has two solutions (—1, —1, —1, —1, —1, 0), (—3, —1, —1, —1, 1, 0) un- 
der (19) and one solution (— 1, — 1, — 1, 1, — 3, 0) under (20). Therefore, there 
are three nonisomorphic (A4(20, 4, 2, 2)'s. For A = 7, d = 2, n = 28 and m = 4, 
there are seven nonisomorphic 0^4(28, 4, 2, 2) 's given by the solutions 

= (-1,-1, -1,-1, 1,1), (-3, -3, -1,-1, 1,0), 
(-5, -1,-1, -1,1,0), (-3, -1,-1, -1,-1,0), 
(-1,-1, -1,-1, -3,0), (-3,-1,-1,1, -3,0), 
(-1,-1,-1,1,-5,0), 



ENUMERATING ORTHOGONAL ARRAYS 



13 



where the first four are solutions under (19) and the remaining three are 
solutions under (20). Similarly, we have the following conclusion for orthog- 
onal arrays of strength 4. There exist no 6X4(48,6,2,4), one nonisomorphic 
0^4(80,6,2,4), three nonisomorphic 6X4(112, 6,2,4) 's and seven nonisomor- 
phic 6X4(144, 6, 2, 4) 's. In fact, the following general result holds. 

Corollary 2. Let g(X,d) be the number of nonisomorphic OA(X2 d ,d + 

2. 2, d)'s with even d and odd X. Then g(X, d) = for A < d — 1 and g(X, d) = 

1. 3, 7 for A = d + 1, d + 3, d + 5, respectively. 

For odd d and odd A, a conclusion similar to Corollary 2 also holds; 
see Corollary 4 in Section 3.3. Nonexistence of an OA(X2 d , d + 2, 2, d) for 
A < d — 1 has previously been established by Blum, Schatz and Seiden [1]; 
see also Theorem 2.29 in [13]. Corollary 2 can easily be verified directly using 
Theorem 1 and can also be obtained from the following Lemmas 7 and 8. 

The complete set of solutions to (18) under either (19) or (20) can, in 
fact, be obtained explicitly in a systematic fashion. For ease of presentation, 
we first introduce some notation. Let Z[a,b] denote the set of integers x 
such that a<x<b and 0[a, b] denote the set of odd integers x such that 
a < x <b. Let S± be the set of solutions to (18) under (19) and S2 be the 
set of solutions to (18) under (20). We then have the following results. 

Lemma 7. For even d and odd X, if A < d — 1, then equation (18) has 
no solution under (19). If A> d+1, then equation (18) has at least one 
solution under (19) and the complete set S\ of solutions is given by 

k G Z[0,(A-d-l)/4], 

u m+l e 0[-(A - 4k) /{m + 1), (A - 4k) /(m - 1)], 

u m G 0[-(X-4k + u m+1 )/m,-\u m+ i\], 

Uj G 0[-(X-4k + u j+ i-\ \-u m+1 )/j,u j+1 ] 

for j = m — 1, m — 2, . . . , 2, 

lil = — (A — 4fc + «2 H 1" Um+l)- 

Lemma 8. For even d and odd X, if X<d + 1, then equation (18) has 
no solution under (20). If X > d + 3, then equation (18) has at least one 
solution under (20) and the complete set S2 of solutions is given by 

k G Z[0,(A-d-3)/4], 

u m G 0[-(A -4k- 2)/{m + 1), (A - 4k - 2)/(m - 1)], 

u m+1 G 0[(m - l)\u m \ — u m — (A — 4k),-\u m \ - 2], 
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Table 1 

The number f(n) of nonisomorphic OA(n, 4, 2, 2) 's when n = 4X is not a multiple of 8. 

Note that f(n) = g(X,2) 



II 


12 


20 


28 


36 


44 


52 


60 


68 


76 


f(n) 


1 


3 


7 


15 


28 


48 


79 


123 


184 


n 


84 


92 


100 


108 


116 


124 


132 


140 


148 


f(n) 


268 


379 


523 


709 


943 


1234 


1594 


2032 


2560 


n 


156 


164 


172 


180 


188 


196 


204 






fin) 


3194 


3946 


4832 


5872 


7082 


8482 


10097 







u m -i G 0[-(\-4k + u m + u m+1 )/(m-l),-\u m \], 

uj G 0[-(A — Ak + Uj + i H \-u m+ i)/j,u j+1 ] 

for j = m — 2, m — 3, . . . , 2, 

ui = -(\-4k + u 2 H \-u m+ i). 

The thorough proofs of Lemmas 7 and 8 involve some tedious, yet rather 
straightforward, algebra and are thus omitted here. These two lemmas, in as- 
sociation with Theorem 1, allow the complete set of nonisomorphic OA (X2 d , 
d + 2, 2,d) , s with even d and odd A to be explicitly and systematically 
obtained through their J*-vectors given by Jt 3 = 2 d uj for j = 1, . . . ,m + 1, 
where (u±, . . . ,u m +i, k) £ S = S\ U S2, with Si and S2 given in Lemmas 
7 and 8, respectively. 

Before moving on, we present in Table 1 an application of Theorem 1 
and Lemmas 7 and 8 to orthogonal arrays of strength 2. Due to spatial 
limitations, we only present the number of nonisomorphic orthogonal arrays 
for 12 < n = 4A < 204 with A = 3, 5, 7, . . . , 51. 

Some results on orthogonal arrays of strength 4 are presented in Ta- 
ble 2. The two cases involving 04(80, 6, 2, 4) 's and 04(112,6, 2,4) 's have 
also been solved in a recent paper by Bulutoglu and Margot [2]. Their re- 
sults on 04(144, m, 2, 4) 's are that there are precisely 20 nonisomorphic 
04(144, 8, 2, 4) 's and there exists no 04(144,9,2,4). 



Table 2 

The number f(n) of nonisomorphic OA(n, 6, 2, 4)'s when n= 16A is not a multiple of 32. 

Note that f(n) = p(A,4) 



/f 


80 


112 


144 


176 


208 


240 


272 


304 


336 


fin) 


1 


3 


7 


14 


26 


46 


77 


123 


190 


n 


368 


400 


432 


464 


496 


528 


560 


592 


624 


fin) 


285 


418 


599 


842 


1163 


1582 


2123 


2813 


3684 
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3.2. Enumerating OA(X2 , d + 2, 2, d)'s for even d and even X. Let A = 
2A*. Since A is even, A* is an integer. Now, consider the equation 

(22) A* + «i + • • • + u m+1 = 2k, 

where k > is an integer and the u/s are integers with \uj\ < X* for all 
j = 1, . . . ,m + 1. Further, the u,-'s must satisfy either (23) or (24): 

(23) tti < • • • < u m < -\u m+ i\, 

(24) Ut < ■■■ <Um-i < -\llm\, U m+1 < -\u m \ - 1. 

Theorem 2. Suppose that d is even and A = 2 A* is also even. We then 
have that: 

(i) each solution (u\, ■ ■ ■ , u m+ \, k) to equation (22) under either (23) or 
(24) determines an OA(X2 d , d + 2, 2, d) with J* -vector given by Jt i =2 d+1 uj 
for j = l,...,m + l; 

(ii) the complete set of nonisomorphic OA(X2 d , d + 2, 2, <i)'s for even d 
and even X is given by collecting the arrays obtained in (\) over all the 
solutions to equation (22). 

We note that, unlike Theorem 1, J tj and Uj are now related through 
Jt j = 2 d+l Uj. The proof of Theorem 2 is quite similar to that of Theorem 1. 
Here we only give an outline of the proof. To prove part (i) of Theorem 2, 
for Jf. = 2 d+1 u-i we need to show that 

/ m+1 \ / m+1 \ 

N s = 2- m [n+Y^h stj J tj = [X*+^h st . Uj \/2 = k- ^ Uj 

\ j=l J \ j=l J j:h stj =-l 

is a nonnegative integer for all s C Z m , which must be true because of 
Lemma 6 and (23) or (24). Part (ii) of Theorem 2 follows from Lemma 3(h) 
and Proposition 1. 

Example 2. Consider 04(8A*,4,2,2)'s. For A* = 1, equation (22) has 
one solution given by (u±, . . . ,u$, k) = (—1,0,0,0,0,0) under (23) and an- 
other solution (0, 0, 0, 0, —1, 0) under (24). For A* = 2, equation (22) has three 
solutions (0,0,0,0,0,1), (-1,-1,0,0,0,0), (-2,0,0,0,0,0) under (23) and 
two solutions (-1,0,0,0,-1,0), (0,0,0,0,-2,0) under (24). When A* = 3, 
equation (22) has the following ten solutions under either (23) or (24): 

(-1,0,0, 0,0,1), (0,0,0, 0,-1,1), (-1,-1, -1,0,0,0), (-1,-1,0, 0,-1,0), 
(-2, -1,0, 0,0,0), (-2,0, 0,0, -1,0), (-1,0, 0,0, -2,0), (-3, 0,0,0, 0,0), 



(0,0, 0,0, -3,0), (-1,-1, -1,-1, 1,0). 
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Similarly, the number of nonisomorphic OA(32X* , 6, 2, 4)'s is 2 and 5 for A* = 
1 and 2, respectively. For A* = 3, the number of nonisomorphic OA(32X* , 6, 2, 4) 's 
is nine, with the corresponding solutions to (22) given by 

(-1,0, 0,0, 0,0, 0,1), (0,0, 0,0, 0,0, -1,1), (-1,-1, -1,0, 0,0, 0,0), 
(-1,-1,0,0, 0,0, -1,0), (-2, -1,0, 0,0, 0,0,0), (-2, 0,0, 0,0,0, -1,0), 

(-1,0, 0,0, 0,0, -2,0), (-3, 0,0, 0,0, 0,0,0), (0,0, 0,0, 0,0, -3,0). 

Except for the two extra zeros in the fifth and sixth positions, these nine 
solutions match the first nine solutions for OA(8X* ,4, 2, 2)'s for A* = 3. These 
results generalize to the situation of any d > 4. 

Corollary 3. Let g(X,d) be the number of nonisomorphic OA(X2 d , 
d + 2, 2,d)'s with even d and even A. Then g(X,d) > 1 for all even X. When 
d = 2, we have g(X, d) = 2, 5, 10 for A = 2, 4, 6, respectively. When d>4, we 
have g(X, d) =2,5,9 for A = 2, 4, 6, respectively. 

One can verify Corollary 3 directly from Theorem 2, or derive it using 
the following Lemmas 9 and 10, which provide all solutions to equation (22) 
and are the respective counterparts of Lemmas 7 and 8 in Section 3.1. Let 
S3 denote the set of solutions to (22) under (23), and S4 the set of solutions 
to (22) under (24). 

Lemma 9. For even d and even X, the complete set S3 of solutions to 
(22) under (23) is given by 

k e Z[0,\*/2], 

u m+1 G Z[-(X* - 2k) /(m + 1), (A* - 2k) /(m - 1)], 

u m £ Z[-(X* -2k + u m+1 )/m,-\u m+1 \], 

Uj £ Z[-(X* - 2k + u j+ i H \-u m +i)/j,u j+ i] 

for j = m — 1, m — 2, . . . , 2, 

ui = -(A* - 2k + u 2 H h u m+ i). 

Lemma 10. For even d and even X, the complete set S4 of solutions to 
(22) under (24) is given by 

k e Z[0, (A* - l)/2], 

u m G Z[-(X* -2k- l)/(m + 1), (A* -2k- l)/(m - 1)], 

u m+1 £ Z[(m - l)\u m \ - Um - (A* - 2k),-\u m \ - 1], 
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Um-i G Z[-(X* -2k + u m + u m+1 )/(m- l),-|n m |], 

Uj £ Z[—(\* -2k + Uj + i H \-u m+1 )/j,u j+1 ] 

for j = m — 2, m — 3, . . . , 2, 

ui = — (A* — 2fc + U2 H 1- n m+ i). 

Lemmas 9 and 10, together with Theorem 2, allow all nonisomorphic 
OA(X2 d ,d + 2,2,d)'s with even d and even A to be explicitly and sys- 
tematically constructed through their J* -vectors given by J t . = 2 d+1 Uj for 
j = 1, ... 1, where (u±, . . . , Um+i, k) £ S3U S4 with £3 and S4 given in 
Lemmas 9 and 10, respectively. 

Table 3 presents the number of nonisomorphic OA(n, 4, 2, 2)'s for 8 < 
n = 8A* < 200. For each n, some of the arrays in Table 3 have strength 3. 
According to a result of Seiden and Zemach ([17]; see also [13], page 35), 
the number of nonisomorphic OA (n, 4, 2, 3) 's is precisely [A*/2] + 1, where 
n = 8A*. If n is a multiple of 16, then exactly one of these strength 3 arrays 
also has strength 4. For 8 < n < 100, Li, Deng and Tang [15] identified the 
nonisomorphic OA(n, 4, 2, 2)'s with distinct confounding frequency vectors 
(CFVs). Comparing their Table 1 with our Tables 1 and 3, we see that the 
CFV does a superb job for the case when n is not a multiple of 8, as it is 
capable of discriminating all the OA(n, 4, 2, 2)'s. However, the CFV is not 
so discriminating when n is a multiple of 8. For example, there are two 
nonisomorphic (9^4(40, 4, 2, 2)'s that have the same CFV, these two arrays 
being given by the solutions to (22), 

(til,..., «5, k) = (-1,-1, -1,-1, -1,0), (-1,-1, -1,-1, 1,1). 

Some results on orthogonal arrays of strength 4 are given in Table 4. Bulu- 
toglu and Margot [2] found that the number of nonisomorphic OA(96, m, 2, 4) 's 
is 4 and for m = 7 and 8, respectively. 



Table 3 

The number f(n) of nonisomorphic OA(n, 4, 2, 2)'s when n = 4X is a multiple of 8. Note 

that f(n)=g(X,2) 



n 


8 


16 


24 


32 


40 


48 


56 


64 


72 


f(n) 


2 


r> 


10 


10 


32 


54 


84 


128 


188 


n 


80 


88 


96 


104 


112 


120 


128 


136 


144 


fin) 


270 


376 


517 


694 


919 


1198 


1543 


1960 


2468 


n 


152 


160 


168 


176 


184 


192 


200 






f(n) 


3072 


3792 


4640 


5636 


6792 


8137 


9682 
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3.3. Enumerating 0A(\2 d , d + 2, 2, d)'s /or odd d and odd A. Consider 
the equation 

(25) A + u\ H h u m+ i = 4k, 

where k and the ttj's are unknown integers, k > 0, tti,. .. , w m are odd with 
\uj \ < A — 2 and it m +i is even with |it m +i| < A — 1. Further, the Uj's must 
satisfy 

(26) ui < ■■■ <u m -i <-\um\, u m+1 <0. 

Theorem 3. Consider OA(\2 d ,d + 2,2,dy s for odd d and odd A. We 
then have that: 

(i) each solution (u\, . . . , u m+ i, k) to equation (25) under (26) deter- 
mines an OA(X2 d ,d + 2,2, d) with J* -vector given by Jt j = 2 d Uj for j = 
1, . . . ,m + 1; 

(ii) the complete set of nonisomorphic OA(X2 d ,d+ 2, 2, c£)'s for odd d and 
odd A is given by collecting the arrays obtained in (\) over all the solutions 
to equation (25) under (26). 

We omit the proof of Theorem 3, which is similar to that of Theorem 1 
and involves the use of Lemma 1, Lemma 3(iv), Lemma 4, Lemma 6 and 
Proposition 2. 

Example 3. Consider CL4(8A, 5,2,3) 's. For A = 1, equation (25) has 
no solution under (26). For A = 3, equation (25) has a unique solution 
under (26) given by (u\, . . . , uq, k) = (— 1, — 1, — 1, — 1, 1, 0, 0). When A = 5, 
there are three solutions (— 3, — 1, — 1, — 1, 1, 0, 0), (— 1, — 1, — 1, — 1, — 1, 0, 0) 
and (— 1, — 1, — 1, — 1, 1, — 2, 0). When A = 7, there are seven solutions: 

(-1,-1,-1, -1,1,0,1), (-5, -1,-1, -1,1, 0,0), (-1,-1, -1,-1, 1,-4,0), 
(-3, -3, -1, -1, 1,0, 0), (-3, -1, -1, -1, -1,0, 0), 

(-3, -1,-1, -1,1, -2,0), (-1,-1, -1,-1, -1,-2,0). 

These generalize to OA(\2 d , d + 2,2, d)'s with A = d - 2, d, d + 2, d + 4 for 
any d > 3. 



Table 4 

The number f(n) of nonisomorphic OA(n, 6, 2, 4)'s when n= 16A is a multiple of 32. 

Note that f(n) = gr(A,4) 



n 


32 


64 


96 


128 


160 


192 


224 


256 


288 


f(n) 


2 


5 


9 


17 


29 


49 


77 


120 


179 


n 


320 


352 


384 


416 


448 


480 


512 


544 


576 


fin) 


265 


380 


539 


747 


1025 


1383 


1848 


2435 


3181 
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Table 5 



The 


number fin) 


of nonisomorphic OA(n, 5, 2, 3)'s 
Note that f(n) = g 


when n = 8A 
(A, 3) 


is not a 


multiple 


of 16. 


n 


24 


40 


56 


72 


88 


104 


120 


136 


152 


fin) 


1 


3 


7 


15 


28 


49 


82 


130 


199 


n 


168 


184 


200 


216 


232 


248 


264 


280 


296 


fin) 


296 


428 


605 


839 


1142 


1530 


2022 


2637 


3399 


n 


312 


328 


344 


360 


376 


392 


408 






fin) 


4336 


5476 


6854 


8509 


10481 


12818 


15573 







Corollary 4. Let g(\,d) be the number of nonisomorphic OA(X2 d , 
d + 2, 2,d)'s with odd d and odd X. Then g(X,d) = for A < d — 2 and 
g(X, d) = 1,3,7 for A = d, d + 2, d + 4, respectively. 

Corollary 4 can be verified directly using Theorem 3, or derived from the 
following Lemma 11, which gives the complete set of solutions to (25) under 

(26) . To present Lemma 11, we need to introduce the notation E[a, b], which 
represents the set of even integers x satisfying a<x<b. 

Lemma 11. For odd d and odd A, if X < d — 2, then equation (25) has no 
solution under (26). If A > d, then equation (25) has at least one solution 
under (26) and the complete set of solutions is given by 

k e Z[0,(A-d)/4], 
u m+1 e E[-{X-d-4k),0], 

u m £ 0[-(X - Ak + u m+ i)/m, (X-4k + u m+ i)/ (m - 2)], 
u m -i £ 0[-(X- 4k + u m +u m+1 )/(m-l),-\u m \], 

Uj G 0[-(X- 4k + u j+1 H \-u m+ i)/j,u j+1 ] 

for j = m — 2, m — 2, . . . , 2, 

ui = -(A - Ak + u 2 H h u m+ i). 

Based on Theorem 3 and Lemma 11, all nonisomorphic OA(X2 d ,d + 
2, 2,d)'s for given odd d and odd A can be obtained. In Table 5, we present 
the number of nonisomorphic (A4(8A, 5, 2, 3)'s for A = 3, 5, ... , 51. 

3.4. Enumerating OA(X2 d , d + 2, 2, (i)'s for odd d and even X. Let A = 
2A*. Since A is even, A* is an integer. Now consider the equation 

(27) A* + ui + • • • + u m+1 = 2k, 
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where k and the Uj's are unknown integers, k > and \uj \ < A* for j = 1, . . . , 
m + 1. Further, the ttj's must satisfy 

(28) < ••■ <it m -i < — \u m \, u m+ i<0. 

Theorem 4. Consider OA(X2 d ,d + 2,2, d)'s for odd d and even A = 2A*. 
We then have that: 

(i) each solution (u\, . . . , u m+ i, k) to equation (27) under (28) deter- 
mines an OA(X2 d , d + 2, 2, d) with J* -vector given by Jtj = 2 d+1 Uj for j = 
1, . . . ,m + 1; 

(ii) the complete set of nonisomorphic OA(\2 d ,d+ 2, 2, cf)'s for oddd and 
even A is given by collecting the arrays obtained in (\) over all the solutions 
to equation (27) under (28). 

Theorem 4 is similar to Theorem 2 in that Jt- and Uj are related through 
Jt j = 2 d+1 Uj. Its proof is also straightforward. Similar to Corollary 3, the 
following result holds. 

Corollary 5. Let g(X,d) be the number of nonisomorphic OA(X2 d , 
d + 2, 2,dys with odd d and even A. Then g(\,d) > 1 for all even A. When 
d = 3, we have g(X, d) = 2, 5, 10 for A = 2, 4, 6, respectively. When d > 5 we 
have g(X, d) =2,5,9 for A = 2, 4, 6, respectively. 

All nonisomorphic CX4(A2 d , d + 2, 2, <i)'s for odd d and even A = 2A* can 
be obtained from Theorem 4 and the following Lemma 12, which gives all 
solutions to (27) under (28). 

Lemma 12. For odd d and even A, the complete set of solutions to (27) 
under (28) is given by 

k G Z[0,X*/2], 

u m+1 G Z[-(A*-2fc),0], 

u m G Z[-(X* - 2k + u m+1 )/m, (A* - 2k + u m+1 )/(m - 2)], 

u m ^i G Z[-(X* - 2k + u m +u m+ i)/(m- l),-|u m |], 

Uj G Z[-(X* -2k + uj + i H \-u m+ i)/j,Uj +1 ] 

for j = m — 2, m — 2, . . . , 2, 

ui = -(A* - 2k + u 2 H h u m+ i). 

Applying Theorem 4 and Lemma 12 to the case d = 3, we have obtained 
all CX4(16A* , 5, 2, 3)'s with A* = 1, . . . , 25. In Table 6, we present the number 
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of nonisomorphic OA(16X* ,5, 2, 3)'s for these values of A*. Some of the arrays 
in Table 6 have strength 4 and the number of such nonisomorphic arrays is 
[A*/2] + 1 ([13], page 35). When n is a multiple of 32, exactly one of these 
strength 4 arrays has strength 5. 

4. Discussion. Using the theory of J-characteristics, we have completely 
solved the problem of enumerating two-level orthogonal arrays of strength 
d with d + 2 constraints. This further demonstrates that J-characteristics 
provide a powerful tool for studying two-level orthogonal arrays. Previously, 
J-characteristics have been successfully used in defining generalizations of 
minimum aberration and in providing statistical justifications for the criteria 
of generalized aberration. 

It is tempting to consider the problem of enumerating two-level orthog- 
onal arrays of strength d with m > d + 3 constraints using the approach in 
this paper. This seems quite nontrivial. Although the general ideas continue 
to apply, useful results like Theorems 1-4 may not exist after all. The fol- 
lowing discussion should shed some light on the issue. First, let us consider 
enumerating orthogonal arrays of strength d with m = d + 1 constraints. 
Based on Lemma 2, we only need to consider the single J-characteristic for 
all m = d + 1 columns. It is also easy to see that we can assume that this 
J value is nonpositive. Using Lemmas 1 and 4, we see that the complete 
set of nonisomorphic OA(X2 d , d + l,2,d)'s is given by the set of solutions 
to the equation A + u = 2k, where k > is an integer and —A < u < 0. 
This provides a very simple derivation for the result of [17] on nonisomor- 
phic OA(X2 d ,d+l,2,dy S . As much as we like this simple derivation, it also 
points to the difficulty arising from trying to solve the problem with d + 3 
constraints. Judging from the huge jump, in terms of difficulty, from the case 
of d + 1 to the case of d + 2, one can only imagine what it will take to solve 
the case of d + 3 completely. In our enumeration of OA(X2 d , d + 2, 2, d)'s, 
Lemma 6 plays an instrumental role. It achieves two goals, one being that it 
substantially narrows down the possible configurations of J-characteristics 
and the other being that we only need to check that N s is a nonnegative 



Table 6 

The number f(n) of nonisomorphic OA(n, 5, 2, 3) 's when n = 8A is a multiple of 16. 

Note that f(n) = g(X, 3) 



n 


16 


32 


48 


64 


80 


96 


112 


128 


144 


f(n) 


2 


5 


10 


19 


33 


56 


89 


138 


207 


n 


160 


176 


192 


208 


224 


240 


256 


272 


288 


f(n) 


303 


432 


606 


832 


1126 


1501 


1975 


2566 


3300 


n 


304 


320 


336 


352 


368 


384 


400 






/(») 


4198 


5293 


6615 


8202 


10092 


12335 


14975 
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integer for s = <j). Although one can derive a result like Lemma 6 for the case 
of m = d + 3, it will not be as powerful as Lemma 6 in achieving the two 
aforementioned goals. Having stated all of the negatives, we recall that ev- 
erything has its positive side. Completely solving the enumeration problem 
for the case of m = d + 3 may just be too ambitious after all. It is, however, 
possible to find a partial solution and this is what we will focus on in our 
future research on this problem. 

Acknowledgments. The authors wish to thank an Associate Editor and 
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